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RESUME 

On  a etudiS  thioriquement  1 'application  du  traite- 
ment  de  signaux  numeriques  ^ un  syst&me  d'alarme  d' intrusion 
et  compare  les  performances  de  quatre  types  differents  de 
f iltres-comparateurs.  Des  simulations  en  ordinateur  nous 
ont  permis  d'evaluer  la  probabilite  de  detection  et  le  taux 
de  fausses  alarmes.  (NC) 


ABSTRACT 

application  of  digital  signal  processing  to 
an  intrusion  alarm  device  is  studied  theoretically.  The  per- 
formance of  four  different  types  of  filter-comparator  units 
is  evaluated  and  compared.  Computer  simulations  are  also 
carried  out  to  obtain  the  probability  of  target  detection  and 
the  false  alarm  rate.  (U) 
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1.0  INTRODUCTION 


This  report  presents  the  results  of  a theoretical  inves- 
tigation of  the  application  of  digital  signal  processing  to  an 
intrusion  alarm  device.  The  investigation  was  carried  out 
in  parallel  with  the  construction  of  a working  model  of  the  in- 
trusion alarm  at  Laval  University  under  a DRB  contract  [1].  The 
theoretical  studies  and  the  experimental  work  involved  in  the 
construction  of  the  alarm  were  both  guided  by  the  findings 
of  an  earlier  DREV  project.  A description  of  the  main  results 
of  this  project  can  be  found  in  a recent  patent  application 
by  Giroux  [2]. 

Although  the  prototype  of  the  intrusion  alarm  operates  in 
the  infrared  spectral  region,  the  concept  described  in  Ref.  2 is  of 
a general  nature.  There,  the  basic  principle  of  the  apparatus  was 
very  clearly  stated:  "In  accordance  with  the  present  invention, 
apparatus  for  detection  of  variations  in  the  radiation  pattern  from 
a defined  space  is  provided  which  comprises  means  for  deriving  a 
digital  signal  from  the  scanned  scene  and  means  for  storing  the 
digital  signal  as  a series  of  discrete  digital  values  in  a digital 
memory.  A digital  comparator  is  provided  for  comparing  the  digital 
output  for  one  scan  with  the  stored  digital  values  in  the  memory 
for  a preceding  scan  of  the  scene.  If  there  is  a discrepancy 
exceeding  some  pre-determined  threshold  value  (to  account  for  noise 
signals,  slow  changes  in  ambient  light  etc)  between  the  digital  values 
for  a given  portion  of  the  scanned  scene  and  the  stored  digital  values 
for  the  same  portion  of  the  scanned  scene  obtained  on  a previous  scan, 
then  an  alarm  is  actuated". 

In  this  report,  we  examine  theoretically  the  application  of 
the  above  principle  to  the  infrared  intrusion  alarm  (IRIA).  To  do  so 
the  first  task  is  to  formulate  the  problem  taking  into  account  the 
practical  requirements  and  constraints.  The  function  of  the  digital 
filter  and  comparator  is  discussed  in  Sec.  2.0,  after  a brief 
description  of  the  general  layout  of  the  IRIA.  Details  concerning 
the  mechanical  and  optical  components  of  the  device  are  given  in 
Ref.  1.  In  Sec.  3.0  the  physical  quantities  based  on  which  various 
possible  filter-comparator  units  are  evaluated  are  studied.  These 
include  the  time  of  response,  the  system  transfer  function,  the 
probability  of  target  detection  and  the  rate  of  false  alarm.  Compa- 
risons between  different  filter-comparator  units  are  made.  To  verify 
the  analytical  results  of  Sec.  3.0  filter-comparator  units  have  been 
simulated  on  a computer.  The  results  of  the  simulation,  the  performance 
characteristics  and  the  relative  merits  of  the  units  considered  are 
described  in  Sec.  4.0.  Section  5.0  summarizes  the  main  results  and 
makes  specific  recommendations  for  the  implementation  of  the  digital 
filter-comparator  unit  to  be  incorporated  to  the  IRIA.  This  work  was 
performed  at  DREV  between  April  1974  and  April  1975  under  PCN  33A07 
(formerly  PCN  15B25,  Project  97-82-06)  "Intrusion  Detection  and  Alarm." 
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2.0  FORMULATION  OF  THE  PROBLEM 


To  define  the  problem  of  signal  processing  associated  with 
the  IRIA,  it  is  useful  to  give  a brief  description  of  the  device.  In 
the  IRIA,  the  incoming  signal  is  obtained  by  scanning  a scene  with  a 
plane  mirror  (Fig.  1).  The  plane  mirror,  driven  by  a motor,  can 
rotate  at  a maximum  angular  speed  of  6.3  rad/s.  The  scene  is 
sampled  at  a rate  of  1000  points  per  revolution.  A pyroelectric  TGS 
infrared  detector  measures  the  infrared  radiation  between  2 and  20  ym 
reflected  by  the  plane  mirror  and  then  focussed  onto  it  by  a concave 
spherical  mirror.  The  instantaneous  field  of  view  of  the  instriiment 
is  approximately  13.2  mrad  by  13.2  mrad.  The  signal  voltage  generated 
at  the  detector  is  converted  to  digital  form  for  each  sampling  posi- 
tion of  the  scanning  mirror. 

For  a given  scan  y (y^O) , a sequence  of  digital  signals 
M (xi,y),  M (X2,  y) , ....  M (x,^,  y)  is  obtained.  The  quantity 
M (x,  y)  is  the  signal  voltage  generated  when  the  intrusion  alarm  scans 
the  spatial  element  with  an  angular  coordinate  x measured  relative 
to  an  arbitrarily  chosen  axis  perpendicular  to  the  axis  of  rota- 
tion of  the  scanning  mirror.  The  value  of  y which  enumerates  the 
scans  is  incremented  by  unity  as  the  IRIA  completes  one  scan.  The 
signal  sequence  M (x,  y)  passes  through  a digital  filter-comparator 
unit  (DFCU)  after  an  initial  period  of  time  necessary  for  stabilizing 
and  establishing  the  background  signal  level.  The  DFCU  decides  on  a 
continuous  basis  whether,  at  any  given  instant,  the  incoming  signal 
M (x,  y)  indicates  a sudden  significant  change  in  the  infrared  radiation 
level  incident  on  the  IRIA  from  the  direction  x.  To  carry  out  this 
task,  information  on  the  background  radiation  and  the  noise  associated 
with  various  sources  of  fluctuations  must  be  collected,  processed, 
stored  and  updated  as  the  scanning  proceeds.  In  other  words,  the 
DFCU  must  be  designed  in  such  a way  that,  at  any  point  (x,y)  in  the 
two  dimensional  lattice  shown  in  Fig.  2,  it  can  decide  whether  M (x,y) 
signals  a significant  change  in  the  infrared  radiation  which  may  be 
attributed  to  an  intrusion  based  on  the  information  obtained  before 
the  arrival  of  M (x,y). 

From  a practical  point  of  view,  it  is  convenient  to  divide 
the  signal  sequence  for  all  angular  positions  up  to  a given  scan 
into  sub-sequences.  Each  sub-sequence  corresponds  to  a given  direc- 
tion X so  that  the  sub-sequences  for  different  directions  are  processed 
in  parallel  independent  of  one  another.  In  situations  where 
the  background  scene  being  scanned  can  be  divided  into  sectors 
within  each  of  which  the  background  and  the  noise  characteristics 
are  the  same  the  amount  of  parallel  processing  can  be  reduced.  For 
our  discussion  we  consider  the  more  general  situation  so  that  the 
DFCU  can  cope  with  changes  in  the  background  in  each  direction  without 
modification.  Consequently,  for  any  given  direction  x,  updated  in- 
formation for  intrusion  detection  must  be  provided  in  a well-defined 
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manner  by  the  DFCU  as  the  latest  incoming  signal  M(x,y)  is  registered. 
Furthermore,  the  required  information  must  be  extracted  from  the  sig- 
nal sequence  M(x,y')  where  y'  = 0,1,2,  »y-l* 

Consider  the  case  when  the  IRIA  is  scanning  a station- 
ary scene,  i.e.,  the  scene  radiance  is  constant  in  time.  If  the 
effect  of  noise  can  be  neglected,  the  digitized  signal  M(x,y)  is 
independent  of  y for  any  value  of  x.  If  no  intrusion  occurs  dur- 
ing the  (y-ljth  scan  the  signal  M(x,y-1)  can  be  used  as  a refer- 
ence level  against  which  M(x,y)  may  be  measured.  In  a real  situ- 
ation, however,  the  noise  from  various  sources  must  be  taken  into 
account.  One  immediate  result  is  that  no  scene  can  be  regarded  as 
strictly  stationary.  Besides,  only  the  variations  in  M(x,y)  not 
attributable  to  the  fluctuations  due  to  noise  are  of  interest. 

To  formulate  the  problem  in  a way  suitable  for  our  subse- 
quent discussions, let  us  assume  that  we  have  N sub-sequences  of 
incoming  signals: 


M(xi,0).  M(xj.l),  M(xj,2), 
M(x2,0),  M(x2.1).  M{x2,2), 


and  M(Xj^,0),  M(Xj^,2) 

where  xj,  X2,...  and  are  the  directions  sampled  in  each  scan. 
Without  intrusion,  each  signal  M(x,y)  may  be  regarded  as  a sum  of 
two  components.  These  are  the  background  signal  and  the  noise 
component  respectively.  Given  a sub-sequence  of  incoming  signals 
.M(x,y')  and  0£y'<y-l  (x  is  fixed  and  y'  varies  between  0 and  y-1), 
and  given  that  no  intrusion  occurs  in  the  corresponding  period  of 
time  (y-l)T  (where  T is  the  time  required  to  complete  one  scan), 
the  problem  is  to  decide  whether  the  signal  M(x,y)  indicates  a 
significant  variation  not  expected  from  the  behaviour  of  the 
signals  in  the  sub-sequence  M(x,y');  0^v'£y-l  which  precedes 
M(x,y).  To  do  so,  the  DFCU  should  be  designed  to  accomplish  two 
things:  (i)  it  should  predict  the  background  signal  component  of 

M(x,y),  from  M(x,y');  O^'^-l;  and  (ii)  it  should  provide  some 
suitable  statistical  measure  of  the  noise  components  in  the  signal 
sequence  so  that  a significance  level  can  be  established.  In 
other  words,  the  problem  is  to  extrapolate  a time  sequence  in  the 
presence  of  noise  and,  at  the  same  time,  measure  the  statistical 
properties  of  the  noise.  To  determine  whether  M(x,y)  indicates  a 
significant  change  in  the  measured  infrared  radiation  level, it  is 
compared  with  the  extrapolated  background  signal.  With  the  known 
statistical  properties  of  the  noise  component, one  can  find  the 
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probability  that  a given  difference  between  the  incoming  signal 
and  the  background  is  caused  by  random  fluctuations. 


3.0  DIGITAL  FILTERING 


The  previous  section  considered  the  problem  of  signal 
processing  in  the  intrusion  alarm.  We  assume  for  our  discussion  that 
the  time  variation  in  the  background  signal  is  slow  so  that  it  can 
be  regarded  as  a constant  over  several  scanning  intervals.  With  this 
assumption,  the  problem  of  extrapolating  the  background  signal 
becomes  essentially  one  of  filtering  out  the  noise  component  from 
the  incoming  signal  sequence. 

In  the  patent  application  [2]  a filter  has  been  suggested 
for  implementation.  In  this  section  we  study  this  filter  (FI,  for 
brevity),  along  with  some  other  possible  alternative  filters,  in  de- 
tail. 

3.1  Criteria  for  filter  comparison 

The  final  objective  for  the  design  of  the  DFCU  is,  of  course, 
the  maximization  of  the  probability  of  detection  of  intrusion  and  the 
minimization  of  the  false  alarm  rate.  At  the  same  time,  subject  to 
the  practical  constraints,  the  IRIA  as  a whole  should  be  simple  and 
inexpensive  to  implement.  To  compare  the  performance  of  different 
filters  it  is  convenient  to  use  the  following  requirements  as  criteria: 

(a)  The  time  of  response  of  the  digital  filter 
should  be  short. 

(b)  The  filter  should  efficiently  filter  out  the 
noise  component  in  the  signal  sequence  M(x,y) 

(c)  The  filtered  signal  sequence  should  follow  the 
slow  variation  (i.e.,  over  time  intervals  >>  T) 
in  the  background  signals  without  appreciable 
distortion. 

We  will  see  later  that  the  filter  which  performs  best  according  to 
these  requirements  also  gives  a higher  detection  probability  and  a 
lower  false  alarm  rate. 

3.2  Choice  of  digital  filters 

In  the  choice  of  filters  we  have  to  restrict  ourselves  to 
those  which  are  simple  to  implement  and  need  a minimum  of  storage 
memory.  Here  we  consider  four  different  types  of  filters:  FI,  F2,  F3 

and  F4.  Some  of  the  properties  of  FI  have  been  given  in  Ref.  2;  the 
impulse  response  function  of  FI  together  with  that  of  the  other  filters 
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are  shown  in  Fig.  3.  If  the  impulse  response  function  of  a given 
filter  is  denoted  by  w (y)  and  the  input  signal  sequence  is 
M(x.y’),  0<y’</,  the  output  of  the  filter  at  time  is 

_ Y 

M(x,y)  = V w(y')  M(x.y-y') 

>^=0 

y 

= H w(y-y')  M(x.y’)  (1) 

y'  = 0 

As  usua^  we  assume  that  the  filter  is  causal,  i.e.  the  output 
signal  M(x,y)  depends  on  the  inputs  from  the  initial  moment  up 
to  the  moment  yT.  To_simplify  our  notations  we  omit  the  vari- 
able X in  M(x,y)  and  M(x,y).  All  the  equations  involving  these 
quantities  are  valid  for  an  arbitrary  but  fixed  value  of  x.  Thus, 
we  replace  Eq.  1 by: 

M(y)  = w(y')  M(y-y')  = ^ w(y-y')  M(y')  (1') 

^=0  y'=0 

The  impulse  response  function  completely  specifies  the  properties 
of  the  filter.  The  analytical  expressions  of  the  impulse  response 
function  of  the  four  types  of  filters  are: 


FI; 

wi(y)  = 

(l-e‘“) 

e (a  > 

0) 

(2) 

F2; 

W2(y)  = 

1/yo 

for  y < 

yo 

= 

0 

otherwise 

(3) 

F3: 

ws(y)  = 

ri  - e 

— 1 e"“y 

for  y < yg 

Ll-e 

■“yol 

= 

0 

otherwise 

(4) 

F4: 

(y)  ^ 

2 

-ay 

= e 

/ 

/ t e 

y'=0 

-“y"  (a  > 0) 

(5) 

Both  FI  and  F4  have  impulse  response  functions  of  infinite 
duration.  While  the  impulse  response  of  FI  decays  exponentially, 
that  of  F4  is  of  the  gaussian  type.  In  Ref.  2,  FI  is  characterized 
by  a constant  K which  is  related  to  a by  the  simple  relation 
K = e°*.  The  impulse  response  functions  of  F2  and  F3  are  of  finite 
duration. 
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Substituting  Eqs.  2 to  5 in  Eq.  l,we  can  write  down 
the  output  of  each  filter  in  terms  of  the  inputs.  However,  it 
is  often  possible  to  write  down  the  output  in  a more  compact 
form  using  the  inputs  as  well  as  the  preceding  outputs.  It  can 
be  readily  verified  that  the  following  recurrence  relations  hold 
for  FI,  F2  and  F3  respectively: 


FI:  M(y)  = (1  - e'“)  M(y)  ♦ e'“  M(y  - 1) 

= (!-!_)  M(y)  * I M(y  - 1) 

K K 


(6) 


F2:  M(y) 


F3:  M(y) 


= M(y)  + TlCy  - 1)  for  y < 

= 7 M(y)  =7  M(y  - y_)  + M(y  - 1)  for  y>y  (7) 

M(y)  + e'“  MCy  - 1)  for  y < yQ 

jwCy)  - e'“^0  M(y  - y^jj  + e‘“  "MCy  - 1) 

for  y>y^j  (8) 


1 - e 


1 - e-^^O 


1 - e 


1 - e-“yo 


Eqs.  6 to  8 provide  the  basic  equations  for  the  implementation 
of  filters  FI,  F2  and  F3.  A similar  equation  does  not  exist  for  F4 
and  it  is  necessary  to  rely  on  Eq.  1 directly.  In  contrast  to  the 
recurrence  equation,  Eq.  1 for  F4  contains  a much  larger  number  of 
terms  and  requires  a correspondingly  larger  memory  for  implementation. 
In  fact,  F4  is  a non-recursive  filter,  while  FI,  F2  and  F3  are 
recursive. 


3.3  Response  time 

The  first  filter  characteristic  we  wish  to  consider  is 
the  response  time.  To  facilitate  the  subsequent  discussion  we 
shall  define  the  response  time  x of  a filter  as  the  time  required 
(in  units  of  the  scanning  time  T)  for  the  filter  output  to  rise 
to  within  2%  of  the  value  of  a time- independent  input  sequence. 

The  need  for  such  a definition  of  x arises  because  its  functional 
form  changes  from  one  filter  to  another.  Since  the  impulse  re- 
sponse functions  of  the  four  types  of  filters  we  consider  are  non- 
negative, the  response  time  is  given  (in  units  of  the  scanning 
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time  T)  by  the  minimum  value  of  y with  which  the  following  inequal- 
ity is  satisfied: 

y 

r w(y')  > 0.98  (9) 

“=0 

For  FI,  the  response  time  can  be  determined  from  the  minimum  value 
of  y so  that: 


^ (1  - e'“)  e'“^'  > 0.98 

y'=0 
or 

> 0.02 

If  T]  is  the  response  time  of  FI  we  obtain: 

. ^ In  0.02 
Tj  = Int  - 

a 


where  Int  [x]  denotes  the  smallest  integer  such  that  Int  [x]  > x. 

From  Eq.  3 it  is  evident,  since  the  impulse  response 
function  of  F2  is  rectangular,  that  the  response  time  T2  is  equal  to 
Yq  foT  Yq4^  50. 

For  F3,the  time  of  response  T3  is  given  by  the  least  value 
of  y for  which  the  following  inequality  is  satisfied 


Hence 


T,  = -1  + Int 


>0.98  (1  - e‘“^0) 

C13) 

r In  (0.98e'“^0  - 0.02)  1 

(14) 

L “ J 

For  a fixed  value  of  T3»more  than  one  combination  of  the  parameters 
a and  y^  may  be  found  so  that  Eq.  14  is  valid. 

The  fact  that  F4  is  non-recursive  makes  it  impossible  to 
write  down  an  exact  expression  for  the  response  time  ti,.  By  defi- 
nition, 14  is  the  minimum  value  of  y which  satisfies  the  inequality: 
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-ay 


i2 


)P=0 


E 

y'=0 


rt2 


“y  > 0.98 


(15) 


To  calculate  T4  approximately,  we  write: 
m 

-an^ 


and 


S(m)  = 

E 

n=0 

m+1 

Si  (m)  H 

E 

n=l 

m 

r 

S (m)  i 

I 

-an^ 


= S(m)  - 1 . 


e dx 


(16) 


(17) 


(18) 


The  integral  S (m)  can  be  expressed  as: 

4/2  am 

S (m)  = ;^  I e 

^2a  *^0 

where  P(x)  is  the  error  integral: 


h _ F 


|P  (-y/^)  - 0.S|  (19) 


P(x)  = £ r dt 


Using  the  approximation: 

S(m)  = i [S(m)  + Sj  (m)] 


and  Eq.  19, we  have: 


S(m)  I P(^2am)  - O.sj  +0.5  [l 

The  above  equation  together  with  the  fact  that : 


(20) 


^-a(m  + 1)2^  ^21) 
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enables  us  to  write  Eq.  15  as: 


^-a(m  + 


1)2 


(22) 


The  least  value  of  m satisfying  Eq.  22  gives  the  time  of  response 
T4.  Equation  22  gives  a sufficiently  accurate  method  for 
calculating  T4  without  having  to  calculate  the  series  in  Eq.  15 
directly  by  term-wise  summation. 

^.4  Frequency  response  and  phase  transfer  functions 

Up  to  this  point,  our  discussion  of  the  various  types 
of  digital  filters  has  been  restricted  to  the  time  domain. 
However,  it  is  instructive  to  study  the  behaviour  of  the  filters 
in  the  frequency  domain.  The  mathematical  tool  used  for  this 
purpose  is  the  Z-transform.  For  any  time  sequence,  iN(y)l  the 
Z-transform  is  defined  by:  ' ' 

QO 

Wn(z)  = ^ N(y).  (23) 

from  which  the  original  sequence  can  be  obtained  by  the  inversion 
formula: 


N(y)  = - /w^(2)  z^'^dz  . (24) 

2iri  J 

unit 

circle 


If  the  impulse  response  function  of  a causal  digital  filter  is  given, 
the  transfer  function  W(z)  of  the  filter  is  defined  as  the  Z-transform: 


W(z)  = 


w(y)  z'^ 


and  the  corresponding  inversion  formula  is: 


w(y) 


f 

W(z)  z^“^dz 


(25) 


(26) 


unit  circle 
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The  properties  of  a given  digital  filter  are  completely  specified 
by  its  transfer  function.  Since  W{z)  is  a complex-valued  function, 
one  may  separate  its  real  and  imaginary  parts: 

W(z)  = (z)  + iW^  (z)  (27) 


R I 

such  that  functions  W (z)  and  W (z)  are  real  valued.  The 
frequency  response  function  R(u))  of  the  filter  is  given  by: 

R(w)  (z  = ^ (z  = e^‘^’^)j  ^ (28) 

where  T is  the  sampling  interval.  The  phase  transfer  function 
(w)  is  given  by  the  formula: 

iuiT. 

' attH- 

From  Eqs.  28  and  29  we  can  see  that  the  frequency  response  and  the 
phase  transfer  functions  are  determined  by  the  behaviour  of  W(z)  on 
unit  circle  of  the  complex  z-plane. 

For  the  filters  FI,  F2  and  F3,it  is  straightforward  to 
calculate  the  transfer  functions  as  the  corresponding  series  (see 
Eq.  25)  can  be  summed  readily.  In  fact,  for  FI: 

00 

Wi(z)  = Sp  Wi(y)  z'^ 


1 ' e ''  1 e z~^ 


(■  • •-)  t (•■•  ■") 


1 - e 


, -a  -1 
1 - e z 


(30) 


For  F2: 


,(z)  = 


%> 


^2  (y)  z~y 
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-h 

y=( 

Similarly, one  can  show  for  F3  that: 
W2(z)  = 


-y  Cl  - z'^0) 

^0  V--') 


(31) 


(1  - 6'°*)  (1  - z'^O) 


(32) 


(1  - e'^^O)  (1  - e-'^z-^) 


-a_-l. 


For  F4  we  are  unable  to  reduce  the  transfer  function  from  the 
definition : 


W4(z)  = S 

y=0 

to  a more  compact  form. 


e ^ z 


-/t 


-ay" 


(33) 


The  frequency  response  functions  for  FI,  F2  and  F3  can 
be  obtained  from  the  transfer  functions  and  Eq.  28.  They  are  given 
by  the  expressions: 


FI:  Ri(w)  = (1  - e'“)y(l  - 2e'“  cosuT  + e ^“) 


F2:  Ra(w)  = (1  = cos  y^  j (1  - cos  wT)  ^ 

(1  - e'“)  (1  - 2e'“^0  cos  y„  toT  + e'^“^0)^  (36) 


(34) 


(35) 


F3:  R3(i^) 


__0 

(1  - e ^^O)  (1  - 20°"  cos  u)T  + e ^“)^ 


The  frequency  response  functions  for  FI,  F2,  F3  and  F4 
are  compared  for  given  times  of  response  in  Figs.  4 to  6. 
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FIGURE  6 - Frequency  response  functions  for  response  time  equal  to  15  T 


4.0  NUMERICAL  SIMULATION  OF  DIGITAL  FILTERS 

To  study  the  efficiency  of  the  filters  in  noise  reduction 
and  their  performance  characteristics  in  terms  of  the  detection 
probability  and  false  alarm  rate,  we  have  carried  out  numerical  simu- 
lations on  a general-purpose  computer.  The  details  of  the  simulation 
are  described  in  a separate  document  [3] ; here,  only  the  main  results 
are  described. 

4.1  Noise  reduction 


As  mentioned  before,  one  of  the  main  objectives  of  intro- 
ducing the  DFCU  is  to  extract  from  a given  incoming  signal  sequence 
the  background  signal  level  for  comparison  with  the  current  data 
point.  The  extent  to  which  a given  filter  achieves  this  can  be 
measured  by  the  corresponding  noise  reduction.  For  this  purpose  we 
may  consider  the  four  types  of  filters  with  response  times  equal 
to  5T,  lOT  and  1ST,  and  two  different  input  signal  sequences.  The 
first  sequence  is  given  by: 

M(y)  = C + n(y)  (37) 

where  C is  a constant  and  n(y)^  a random  gaussian  noise  with 
a standard  deviation  equal  tO'^C.  To  compare  the  noise  in  the  out- 
put sequence  with  that  in  the  input  sequence  the  ratio: 


f 


IJNCI.ASSU-  li;i) 
1<) 


var  M(y) 
var  M(y) 


y'=yo 


^ M(y')  - < M(y')> 

^M(y')  -<M(y')> 


(38) 


can  be  calculated.  The  initial  point  y^.  in  the  sununations  of  Eq.  38 
is  chosen  so  that  it  is  larger  than  the  response  time  and  the  final 
point  N+y„  is  chosen  so  that  the  results  obtained  do  not  show  signif- 
icant variation  when  additional  terms  are  included  into  the  summations. 
The  angular  bracket  < > indicates  the  mean  value  of  the  quantity 

enclosed  taken  in  the  time  interval  under  consideration. 


Another  comparison  can  also  be  carried  out  for  the  four 
types  of  filters  using  two  signal  sequences: 

Mi(y)  = C sin  (Y  = 30) 

and  " 

MaCy)  = C sin  + n(y) 

Y 

With  the  filtered  output  the  mean  square  error,  (MSE)  as  defined 
below,  is  calculatecl: 


( 

[ M2  (y')  - M (y')V 

y’=yo 

^ 1 / 

Mi^  (y) 

(39) 

y'=yo 

The  results  of  these  calculations  are  given  in  Tables  I to  III  which 
show,  that  for  an  equal  response  time,  F2  and  F3  (for  the  best 
combination  of  the  parameters  a and  Y^)  give  about  the  same  noise 
reduction  with  a constant  or  sinusoidal  background.  Although  F4 
does  not  reduce  the  noise  as  effectively  as  F2  or  F3,  it  does 
give  better  results  than  FI.  However,  from  the  phase  transfer 
functions  of  these  filters  (Tables  I to  III)  it  can  be  seen  that, 
because  of  their  larger  phase  shifts,  F2  and  F3  will  tend  to  intro- 
duce more  errors  than  FI  when  a quickly  varying  background  is 
present. 
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Apart  from  separating  the  background  signal  from  noise  it 
is  desirable  that  the  DFCU  should  establish  and  update  the  detection 
threshold.  Since  the  background  is  assumed  to  be  slowly  varying 
we  may  calculate  the  root -mean-square  noise  level  o from  the 
equation 


'^(y) 


t [ 


M(y')  - M(y' 


y'=y-yr 


or  in  an  equivalent  but  recursive  manner; 


o(y)  (y  - 1) 


‘ 7*.  [' 


- M(y  - 1)  C41) 


The  choise  of  a suitable  value  of  N in  Eq.  40  is  determined  by  the 
consideration,  on  the  one  hand,  that  it  should  be  sufficiently  large 
so  that  o(y)  is  not  sensitive  to  the  change  in  the  length  of  the 
time  interval  in  the  summation  and,  on  the  other  hand,  that  it 
should  be  kept  as  low  as  possible  to  reduce  the  'warm-up'  time  of 
the  whole  system. 

To  update  the  value  of  o(y)  we  may  either  calculate  a sepa- 
rate a(y)  from  the  incoming  signal  while  the  o(y)  already  available 
is  kept  until  the  calculation  is  completed.  A simpler  method  is  to 
update  a according  to  the  equation; 


o(y)  = 


J o (y)  . j7 

▼ N+1  N 


M(y)  - M(y  - 1)1 


Both  methods  have  been  tested  and  they  gave  essentially  the  same 
results. 

4,2  Detection  probability  and  false  alarm  rate 


With  a given  value  of  o(y),an  incoming  signal  M(y)  can 
be  compared  to  the  filtered  signal  M(y  - 1)  to  determine  which  of 
the  following  two  relations  is  satisfied; 


|M(y)  - MCy  - 1)1  < ko 

|M(y)  - "MCy  - 1)1  > ko 


(43  a) 
(43  b) 


To  assess  the  performance  of  the  different  types  of  filters 
we  may  assume  a constant  background  with  a gaussian  noise  added  to 
a randomly  occurring  signal  of  such  a magnitude  that  the  signal- 
to-noise  ratio  is  equal  to  three.  The  detection  probabilities 
and  the  false  alarm  rates  were  found  for  different  values  of  k in 
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relations  43  a and  6.  The  results,  which  are  averages  over  eight 
trials,  are  shown  in  Table  IV.  For  all  the  filters^ the  probabilities 
of  detection  and  the  false  alarm  rates  are  high  for  small  values  of  k. 
Increasing  the  value  of  k lowers  the  false  alarm  rate  and  the  detec- 
tion probability  at  the  same  time.  The  differences  in  the  performance 
of  the  various  types  of  filters  can  be  compared  by  calculating  the 
number  of  the  random  signals  detected  minus  the  number  of  false  alarms 
in  a given  period  of  time  (see  Table  IV) . It  should  be  noted  that 
different  values  of  k are  required  to  optimize  the  different  filters. 


TABLE  I 


RELATIVE  PERFORMANCE  OF  FILTERS 
(RESPONSE  TIME  5 T) 


V I Iter 
boscript ion 


Coefficients 


Transfer  function 

A = normalized  amplitude 
P s phase  (degrees) 

Y » Period 


Y* 


20T 


30T 


40T 


SOT 


Noise 

reduction 

VAR  M (y) 

VAR  M (y) 


Mean 

square 

error 

V 


FI 


a a 0.652 


A=  0.905 
H=  -19.57 


0.955 


0.973 


0.983 


-12.45 


-9.53 


-7.69 


0.337 


2.047 


= 5.000 


A=  0.904 
P=  -36.00 


0.957 


0.975 


0.984 


-24.00 


-18.00  -14.40 


0.237 


1.335 


F3 


a = 0.300 


Yg  = 6.000 


A=  0.877 
P=  -29.98 


0.951 


0.968  0.982 


-19.81 


-14.93  -11.97 


0.245 


1.406 


F3 


a = 0.500 


yp  = 7.000 


A=  0.890 
P=  -22.66 


0.955 


0.970  0.983 


-15.51 


-11.75  -9.45 


0.289 


1.712 


F3 


a = 0.600 


Yg  = 8.000 


A=  0.903 
P=  -19.28 


0.957 


0.974  0.984 


-13.38 


-10.17 


-8. 19 


0.320 


1 .927 


F4 


a = 0.150 


A=  0.934 
P=  -20.55 


0.973 


0.983  0.990 


-13.84 


-10.42  -8.35 


0.294 


1.771 


UNCLASSIFIED 

20 


TABLE  III 


Rl.LATUh  PERFORMANCE  OF  FILTERS 
(RESPONSE  TIME  15  T) 
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5.0  CONCLUSIONS 


From  the  noise  reduction  simulation  and  the  calculation  of 
the  detection  probability  and  the  false  alarm  rate  discussed  above, 
we  can  see  that  among  the  four  types  of  filters  F2  gives  the  best 
performance.  The  results  obtained  for  F.5  are  quite  close  to  those 
obtained  for  F2,  and  those  for  F4  show  that  it  is  the  least 
effective.  For  the  purposes  of  building  a testing  model  it  seems 
only  necessary  to  consider  FI  and  F2.  The  great  advantage  of  FI, 
in  spite  of  its  poorer  performance,  is  the  simplicity  of  its  imple- 
mentation. Its  characteristics  are  completely  determined  by  a 
parameter  K,  which  can  easily  be  varied.  For  F2,  on  the  other 
hand,  the  only  variable  parameter  is  the  memory  capacity  once  the 
number  of  sampling  points  per  scan  is  fixed.  This  is  a significant 
drawback  in  comparison  with  FI. 

The  results  presented  in  Sec.  4.0  for  the  detection  probabil- 
ity and  false  alarm  rate  also  show  that  F2  gives  a better  performance 
than  FI.  To  reduce  the  false  alarm  rate  without  decreasing  at  the 
same  time  the  detection  probability,  we  may  decide  to  trigger  an 
alarm  only  after  observing  an  anomalous  variation  in  the  incoming 
signals  over  a number  of  consecutive  scans  in  a correlated  space 
sector.  To  do  this, only  minor  modifications  are  required  in  the 
signal  processing. 

In  this  study  we  have  not  considered  the  question  of  making 
use  of  more  specific  information  on  targets  and  backgrounds.  This 
is  one  area  in  which  further  work  would  help  to  enhance  the  discrim- 
inating power  of  the  intrusion  alarm. 
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